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SUMMARY 

Surface segregation profiles of binary (Cu-Ni, Au-Ni, Cu-Au) and ternary (Cu-Au-Ni) alloys are deter- 
mined via Monte Carlo-Metropolis computer simulations using the BFS method for alloys for the calculation 
of the energetics. The behavior of Cu or Au in Ni is contrasted with their behavior when both are present. The 
interaction between Cu and Au and its effect on the segregation profiles for Cu-Au-Ni alloys is discussed. 


INTRODUCTION 

Select binary systems have had their surface compositions thoroughly investigated both by experimental 
and by theoretical techniques, as reviewed in reference 1. Recent work also includes similar research for 
combinations of three or more elements (refs. 2 to 8), ranging from metal-semiconductor surfaces (ref. 9) to 
alloy impurity cosegregtion (refs. 10 and 11). The theoretical determination of temperature-dependent 
segregation profiles constitutes a severe test for any modeling effort, as it requires a great deal of accuracy 
from the model used in describing the simultaneous interactions of several atomic species. In this work, we 
present an extensive study of binary and ternary systems of Ni, Cu and Au based on large-scale temperature- 
dependent Monte Carlo/Metropolis simulations using the Bozzolo-Ferrante-Smith (BFS) method for alloys 
(refs. 12 to 21) for computing the energy, and compare with experiment in those cases where data are 
available (refs. 22 to 26). 


THE BFS METHOD AND THE MONTE CARLO/METROPOLIS ALGORITHM 

Since its inception a few years ago, the BFS method has been applied to a variety of problems, ranging 
from bulk properties of solid solution fee alloys (ref. 13) and the defect structure in ordered bcc alloys (refs. 14 
and 15) to more specific applications including detailed studies of the structure and composition of alloy sur- 
faces and surface alloys (ref. 16). 

In what follows, we provide a brief description of the operational equations of BFS. The reader is encour- 
aged to seek further details in previous papers where a detailed presentation of the foundation of the method, 
its basis in perturbation theory and a discussion of the approximations made are clearly shown (refs. 12 to 18). 

The BFS method provides a simple algorithm for the calculation of the energy of formation of an arbitrary 
alloy (the difference between the energy of the alloy and that of its individual constituents). In BFS, the 
energy of formation AH is written as the superposition of elemental contributions e i of all the atoms in the alloy 

i i 

where E\ is the energy of atom i in the alloy and E, is the corresponding value in a pure equilibrium mona- 
tomic crystal. In principle, the calculation of AH would simply imply computing the energy of each atom in its 
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equilibrium pure crystal and then its energy in the alloy. In BFS, beyond directly computing the difference e { 
for each atom in the alloy, a two-step approach is introduced for such a calculation in order to identify contri- 
butions to the energy due to structural and compositional effects. Therefore, e t is broken up in two separate 

contributions: a strain energy (ef) and a chemical energy -ef -ef° ) While there is a certain level of 

arbitrariness in how this separation is implemented, it is only meaningful when a good representation of the 
actual process is obtained by properly linking both contributions. This is achieved by recoupling the strain and 
chemical contributions by means of a coupling function, properly defined to provide the correct asymptotic 
behavior of the individual components. Each individual contribution e i can therefore be written as 

e , = 4 + g( e f - 4 ° ) (2 

s 

The BFS strain energy contribution e i is defined as the contribution to the energy of formation from an 
atom in an alloy computed as if all the surrounding atoms were of the same atomic species, while maintaining 
the original structure of the alloy. To visualize this concept, figure 1(a) represents the atom in question (identi- 
fied with an arrow) in an equilibrium position in its ground state crystal (arbitrarily represented by a simple 
cubic lattice). Figure 1(b) shows the same atom in the alloy being studied (also arbitrarily represented by a 
different crystallographic symmetry). Two things can be different between the reference crystal and the alloy. 
First, atoms of other species may occupy neighboring sites in the crystal and, second, the crystal lattice may 
not be equivalent in size or structure to that of the ground state crystal of the reference atom. In figure 1(b), 
the different atomic species are denoted with different symbols from that used for the reference atom, and the 
differences in size and/or structure are denoted with a schematically different atomic distribution as compared 
to the ground state crystal shown in figure 1(a). The BFS strain energy accounts for the change in energy due 
only to the change in geometrical environment of the crystal lattice (from figs. 1(a) to (b)), ignoring the addi- 
tional degree of freedom introduced by the varying atomic species in the alloy. In this context, figure 1(c) 
shows the environment ‘seen' by the reference atom when computing its BFS strain energy contribution. The 
neighboring atoms conserve the sites in the actual alloy (fig. 1(b)), but their chemical identity has changed to 
that of the reference atom (fig. 1(a)) thus simplifying the calculation to that of a single-element crystal. The 
BFS strain energy term represents the change in energy of the reference atom in going from the configuration 
denoted in figures 1(a) to (c). In this sense, the BFS strain energy differs from the commonly defined strain 
energy in that the actual chemical environment is replaced by that of a monatomic crystal. Its calculation is 
then straightforward, even amenable to first-principles techniques. In our work, we use Equivalent Crystal 
Theory (ECT) (refs. 19 and 20) for its computation, due to its proven ability to provide accurate and computa- 
tionally economical answers to most general situations. In all cases considered in this work, a rigorous 
application of ECT is reduced to that of its two leading terms, which describe average density contributions 
and bond-compression anisotropies. We neglect the three- and four-body terms dealing with the bond angle 
and face-diagonal anisotropies. 

The chemical environment of atom i is considered in the computation of ef , the first term in the total 
BFS chemical energy contribution, where the surrounding atoms maintain their identity but are forced to 
occupy equilibrium lattice sites corresponding to the reference atom /. Following the convention introduced 
in figures 1 and 2(a) shows the reference atom in the actual alloy (similar to fig. 1(b)), while figure 2(b) indi- 
cates the atomic distribution used in computing the BFS chemical energy ef (note that the lattice used in 
fig. 2(b) corresponds to that of the ground state crystal of the reference atom, as shown in fig. 1(a)). The total 
BFS chemical energy is then the difference between the energy of the reference atom in figure 2(b), ef, and 
its energy in its ground state crystal (fig. 1(a)). Building on the concepts of ECT, a straightforward approach 
for the calculation of the chemical energy is defined, properly parameterizing the interaction between dissimi- 

lar atoms. The second contribution to the chemical energy e- 0 (fig. 1(a)) is included in order to completely 
free the chemical energy from structural defects, taking into account the possibility that the reference atom 
is not in a full-coordination environment. This is accomplished by recomputing the contribution ef defined 
before, but once again assuming that all atoms are of the same species as the reference atom. 

As mentioned above, the BFS strain and chemical energy contributions take into account different effects, 
i.e., geometry and composition, computing them as isolated effects. A coupling function, g jt restores the rela- 
tionship between the two terms. This factor is defined in such a way as to properly consider the asymptotic 
behavior of the chemical energy, where chemical effects are negligible for large separations between 


i 
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dissimilar atoms. Within the framework of this discussion, the total BFS contribution e, of each atom in the 
alloy can be graphically depicted by the combination of strain and chemical effects shown in figure 3. 

In what follows, we provide the basic operational equations needed to compute each one of the terms 
introduced above. The BFS strain energy contribution ef is obtained by solving the ECT perturbation equation 


NR^e~ a ' Rx +MR?e 


-(a,+l/A/)«2 _ 


-11 


,Pi d 


►%)k/ 


(3) 


where N and M are the number of nearest- and next-nearest neighbors respectively, and where p, I , a and A are 
ECT parameters that describe element / (see ref. 17 for definitions and details), r denotes the distance 
between the reference atom and its neighbors, S(r) describes a screening function (ref. 19) and the sum runs 
over nearest- and next-nearest neighbors. This equation determines the lattice parameter of a perfect equiva- 
lent crystal where the reference atom i has the same energy as it has in the geometrical environment of the 
alloy under study. R , and R, denote the nearest- and next-nearest neighbor distances in this equivalent crystal. 

Once the lattice parameter of the (strain) equivalent crystal, a s , is determined, the BFS strain energy con- 
tribution is computed using the universal binding energy relation of Rose et al. (ref. 20), which contains all the 
relevant information concerning a single-component system: 



(4) 


where /•.[- is the cohesive energy of atom / and where the scaled lattice parameter af is given by 


s* 


= q( a ? -a'e)/li 


(5) 


where q is the ratio between the equilibrium Wigner-Seitz radius and the equilibrium lattice parameter o'. 

The BFS chemical energy is obtained by a similar procedure. As opposed to the strain energy term, the 
surrounding atoms retain their chemical identity, but are forced to be in equilibrium lattice sites of an equilib- 
rium (otherwise monatomic) crystal i. The BFS equation for the chemical energy is given by 

NRf' e~ a ' Ri +MR^e- {ai+l/Xi)R2 = j (6) 


where N jk and A/,, are the number of nearest- and next-nearest neighbors of species k of atom i. The chemical 
environment surrounding atom i is reflected in the parameters a, k , given by 

ik ~ ®i T ^ki 


where the BFS parameters A (a perturbation on the single-element ECT parameter a , ) describe the changes of 
the wave function in the overlap region between atoms i and k. Once equation (6) is solved for the equivalent 
chemical lattice parameter af, the BFS chemical energy is then 



where = 1 if af * > 0 and y = -l if af *< 0. The scaled chemical lattice parameter is given by 

of* =q(af -a'e)/ 1 ! 


( 8 ) 


(9) 
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Finally, as mentioned above, the BFS chemical and strain energy contributions are linked by a coupling 
function g, which describes the influence of the geometrical distribution of the surrounding atoms in relation to 
the chemical effects and is given by 


-a? 

Si = c ' (10) 

c* 

where the scaled lattice parameter a x is defined in equation (5). 

In this work we used the BFS interaction parameters A determined following the procedure outlined in ref- 
erence 13. The pure element parameters a e , £ f , /, a, A and the BFS parameters A PdCu and A CuPd used in this 
study are listed in Table I. 

Large scale simulations were performed, covering a wide range of concentration for all the alloys consid- 
ered, using the Monte Carlo-Metropolis algorithm and the BFS method for the calculation of the energetics. 
The temperature treatment is simulated by means of an algorithm where pairs of atoms of different species are 
allowed to switch places. The switch is accepted if it lowers the energy, or with probability 
e*/>(-AE/kT) otherwise, where k is Boltzmanns constant, T is the temperature and AE is the difference in en- 
ergy of the cell before and after the switch. The changes in atomic distribution are allowed to continue until 
the total energy of the cell stabilizes. 


RESULTS AND DISCUSSION 

The segregation profiles for several Cu-Ni, Cu-Au and Au-Ni alloys were determined. The results for some 
selected concentrations for Cu-Ni and Au-Ni alloys are displayed in figure 1. While the Monte Carlo simul- 
ations include an appropriate temperature treatment for these alloys, there are restrictions in the calculation 
directed to simplify the numerical task to a reasonable level. Perhaps, the most important approximation re- 
lates to the use of a rigid lattice with a fixed lattice parameter throughout the computational cell, thus disal- 
lowing for multilayer or individual relaxations. In most cases, there is no experimental knowledge of the lattice 
parameters for the bulk ternary alloys. We therefore used approximate values based on the BF rule 
(ref. 21). 

For Cu-Ni, it is observed that Cu always segregates to the surface plane, and to the second layer as well, 
as the Cu concentration (x Cu ) increases. The surface concentration of Cu attains a maximum even for small 
values of bulk Cu concentration, in excellent agreement with experiment and other theoretical models (refs. 1 
to 12, and 22). The monotonic profile characteristic of low Cu concentrations is broken for x Cu - 10 percent, 
where small fluctuations in Cu and Ni concentration per plane become apparent, without any differences in the 
composition of the top surface planes. This behavior progresses to a relatively smooth profile, with less pro- 
nounced oscillations for Cu-rich alloys. It should be noted, however, that the apparent oscillations in the pro- 
files shown in figures 4(a) to (d) most likely represent fluctuations in concentration, highly dependent on the 
temperature treatment in the simulation, rather than actual oscillatory patterns. For Au-Ni, Au segregation to 
the surface plane is observed in all cases, in agreement with experiment (refs. 1 and 23), with the additional 
feature of clearly defined oscillations. While the abundance of Au in the surface has been found experi- 
mentally (ref. 22), no other theoretical or experimental evidence exist to support or contradict the behavior 
observed in figures 4(e) to (h) for the first few planes immediately below the surface plane. Cu-Au profiles 
show a small departure from the corresponding ordered bulk structure characteristic of these alloys beyond the 
top surface planes. It is observed that the surface plane is always rich in Au, in agreement with experiment 
(refs. 24 and 25). In particular, the lCu-99Au case shown in figure 5(a) shows evidence of two competing seg- 
regation mechanisms, where Au populates the surface plane and Cu segregates to the plane immediately be- 
low the surface. This result is in agreement with the experimental finding of full occupancy of Au atoms in the 
first layer inducing the depletion of Au atoms in the layer below (refs. 24 to 26). 

The successful comparison of the results for binary systems with available experimental data and, in par- 
ticular, the examination of three substantially different behaviors for each of the three binary alloys studied, 
provides a solid foundation for extending the approach to model three-component systems. With the solid foun- 
dation developed with these examples, we now study the Cu- Au-Ni system and apply the same method-ology 
to determine the surface composition of a variety of Cu-Au-Ni alloys. 

While both Au and Cu show' a tendency to segregate to the outermost layers in Au-Ni and Cu- Ni alloys, 
the interaction between Cu and Au in Cu-Au alloys indicates that the competition between Cu and Au in 
Ni-rich ternary systems, studied in figure 6, might not be anticipated from the behavior of binary systems 
alone. The profile shown in figure 6(a) indicates that, for the same (low) Cu and Au concentration, 
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Cu segregates to the surface, as is the case when the concentration of Cu exceeds that of Au (figs. 6(a), (c), 
and (g)). Not surprisingly, even for small amounts of Au the tendency for Au surface segregation rapidly over- 
comes that of Cu, to the point where there is a direct competition for populating the surface plane between the 
two components in 4Cu-lAu- 95Ni (fig. 6(c)). For Au concentrations greater than 1 percent, the behavior of Au 
is dominant, inducing segregation patterns characteristic of Au-Ni or Au-Cu. As seen in figures 6(d), (e), and 
(g), Cu is lumped in the few layers below the Au-rich top layers. The best example of this behavior is seen in 
figure 6(d), where the first three layers have a population of Au and Ni similar to that seen in the binary case 
(fig. 5(b)), whereas Cu dominates the concentration pattern of the next few layers in a manner similar to that 
seen in figures 4(b) and (c). Another remarkable case corresponds to the 19Cu-lAu-80Ni where, with the ex- 
ception of the top two layers, the concentration profile matches that seen for 20Cu-80Ni in figure 4(d). How- 
ever, the population of the top two layers is a clear demonstration of the competition between the surface 
segregation properties of Cu or Au, yielding a mixed Cu-Au surface plane. In summary, features that are char- 
acteristic of the binary systems (Au-Ni, Cu-Ni) directly compete in the ternary case, with slight modifications 
due to the interactions that arise between Cu and Au in the surface region. 


CONCLUSIONS 

The calculation of segregation profiles for binary and ternary alloys using Monte Carlo techniques and the 
BFS method for the energetics allows for a direct comparison between individual and collective behaviors of 
each one of the constituent elements. Due to the diverse behavior observed in the binary systems, the method- 
logy was applied to the Cu-Au-Ni system for w'hich no experimental data exists. The possibility of studying 
complex systems with such computational ease could provide a useful tool to assist in the development of ex- 
perimental work on multicomponent systems. 
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TABLE I.— EXPERIMENTAL RESULTS USED FOR THE DETERMINATION OF THE 
ECT AND BFS PARAMETERS 4 



Experimental results 


ECT paramters 


Latice param- 
ter, 

A 

Cohesive en- 
ergy, 
eV 

Bulk modulus. 
GPa 

P 

a 

(A -') 

?. 

(A -') 

1 

(A) 

Pd 

3.890 

3.94 

195.77 

8 

3.612 

0.666 

0.237 

Cu 

3.615 

3.50 

142.08 

6 

2.935 

0.765 

0.272 

[ BFS parameters j 

W--Q-QWS 



= -0.0431 



'The last four columns display the resulting equivalent crystal theory (ECT) (ref. 26) parameters: p is related to the 
principal quantum umber n for the atomic species considered (p = 2n - 2), a parameterizes the electron density in the 
overlap region between two neighboring atoms, X is a screening factor for atoms at distances greater than nearest- 
neighbor distance and / is a scaling length needed to fit the lattice parameter dependence of the energy of formation with 
the universal binding energy relationship of Rose et al. (ref. 27). The last row displays the BFS parameters A^ and A cupd 
used in this work. 



Figure 1 . — (a) A pure, equilibrium crystal (reference atom denoted by the arrow), 
(b) a reference atom (denoted by the arrow) in the alloy to be studied (atoms 
of other species denoted with other shading) and (c) the same reference atom 
in a monatomic crystal, with the identical structure of the alloy to be studied, 
but with all the atoms of the same atomic species as the reference atom, for 
the calculation of the strain energy term for the reference atom. The strain 
energy is the differencein energy of the reference atom between (c) and (a). 
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Figure 2 — (a) The reference atom (denoted by the 
arrow) in the actual alloy environment and (b) the 
reference atom surrounded by a chemical environ 
ment equivalent equivalent to that in (a) but with 
the different neighboring atoms occupying equili- 
brium lattice sites corresponding to the ground 
state of the reference atom. 



Figure 3.— Schematic representation of the BFS contributions to the total energy of formation. The left hand 
side represents the reference atom (denoted by an arrow) in an alloy. The different terms on the right hand 
side indicate the strain energy (atoms in their actual positions but of the same atomic species as the 
reference atom), the chemical energy term (atoms in ideal lattice sites) and the reference chemical energy 
(same as before, but with the atoms retaining the original identity of the reference atoms). 
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(a) 1Cu-99Ni 


(b) 5Cu-95Ni 


(c) 10Cu-90Ni 
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(e) 1Au-99Ni 


(d) 20Cu-80Ni 


(f) 5Au*95Ni 


(g) 10Au-90Ni 


(h) 20Au-80Ni 


Figure 4.— Segregation profiles of xCu-yNi (a-d) and xAu-yNi (e-h) for (x,y) = (1,99), (5,95), (10,90), and (20,80). 
The horizontal axis indicates the atomic planes with the first one corresponds to the surface plane. The 
vertical axis denotes the population of each 72-atom plane. The purple solid line corresponds to Ni while 
the green dash line corresponds to Au (e-h) and the red dot-dash line to Cu (a-d). 
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(a) 1Cu-99Au (b) 5Cu-95Au 



Figure 5.— (a-g) Segregation profiles of xCu-yAu for (x,y) = (1,99), (5,95), (25,75), (50,50), 
(75,25), (95,5) and (99,1). The horizontal axis indicates the atomic planes with the first 
one corresponds to the surface plane. The vertical axis denotes the population of each 
72-atom plane. The green dash line corresponds to Au and the red dot-dash line to Cu. 
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(a) 1Cu-1Au-98Ni (b)1Cu-4Au-95Ni 




(c) 4Cu-1Au-95Ni 



(e) 10Cu-10Au-80Ni 





Figure 6 — Segregation profiles of (a-e) xCu-yAu for (x,y) = (1,1,98), (1,4,95), (4,1,95), 
(5,5,90), (10,10,80), (1,19,80), and (19,1,80). The horizontal axis is labeled by the number 
of planes. The vertical axis denotes the population of each 72-atom plane. The purple 
solid line corresponds to Ni, the green dash line to Au, and red dot-dash line to Cu. 
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